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Abstract 

We give the analytical expressions and the numerical values of radiative correc- 
tions to the covariant derivative operator on the static quark line, used for the lattice 
calculation of the Isgur-Wise form factors ri/2(l) and t^/2{1). These corrections in- 
duce an enhancement of renormalized quantities if an hypercubic blocking procedure 
is used for the Wilson line, while there is a reduction without such a procedure. 

PACS: 12.38.Gc (Lattice QCD calculations), 12.39.Hg (Heavy quark effective theory), 13.20. He (Lep- 
tonic/semileptonic decays of bottom mesons). 



1 Introduction 



In a previous paper [1] we proposed a method to compute on the lattice, in the static 
hmit of HQET, the Isgur-Wise form factors Ti/2(1) and T3/2(1) which parameterize de- 
cays of B mesons into orbitally excited (P wave) D** charmed mesons. Keep in mind 
that the zero recoil is the only definite limit of HQET on the lattice, because the Eu- 
clidean effective theory with a non vanishing spatial momentum of the heavy quark is not 
lower bounded. Then it reveals impossible to calculate directly the r^'s from the currents 
because the matrix elements vanish at zero recoil. To compute ri/2(l) and T3/2(1), wc 
proposed to evaluate on the lattice the matrix elements {HQ{v)\h{v)jij5Djh{v)\H{v)) and 
{H2{v)\h{v)jij5Djh{v)\H{v)), using the following equalities: 

{H;{v)\h{v)^a^Djh{v)\H{v)) ^iQij {Mh^ - Mn) rx (1), (1) 
{H;{v)\h{v)^anD,h{v)\H{v)) = -tV3{MH* - Mh) r|(l)e^ , (2) 

where Di is the covariant derivative {Di = di + igAi), Mh, Mh* and Mh* are the mass of 
the 0^, 0^ and 2+ states respectively, and e*j is the polarization tensor. These relations 
are defined between renormalized quantities. Then we have to renormalize the matrix 
element of the derivative operator computed on the lattice. We explained that power 
and logarithmic divergences are not to be feared in the zero recoil limit. However finite 
renormalization is present and wc want to establish the one-loop contributions to the 
derivative operator with the hypercubic blocking [2] of the Wilson line. 

We have to renormalize and to match onto the continuum the bare operator Ofj — 
h^lil^Djh^ , where is the bare heavy quark field. We choose the MOM scheme whose 
renormalization conditions are: 1) the renormalized heavy quark propagator is equal to 
the free one, and 2) the renormalized vertex function taken between renormalized external 
legs is the tree level one. 

The rhs of equations (1) and (2) are independent of the renormalization scale /x. Indeed, 
on the one hand, Mh^ —Mh = Ao+ — Aq- and Mh* —Mh = A2+ — Aq- are physical quantities. 
On the other hand. 



^mB-mDo 

F^" = ^/^iTTTCf (//, w)a^' + V^^I^ [Clin, w)v'' + Clin, w)v'^] , 

where \^w^ — la^ = (v — v')^, and the Cf's are the matching coefficients between the QCD 
operator cj^j^b and the HQET operators Cy/j^j^by, c^'V^j^by and c^/v'^j^by respectively. 



-n/2{f^,l)V2C!{f,,l)a^. (3) 
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Cfifi, 1) = Cf{l) [3] (Cf (1) =r]A = 0.986 ± 0.005 [4]) and the Ihs of (3) is also indepen- 
dent of ii: thus Ti/2{fi, 1) = ri/2(l). We can use the same argument to prove the scale 
independence of T3/2(/i, 1). Consequently the scale /i will be omitted in the following. 

It is well known that the heavy quark self-energy diverges linearly in 1/a [5], so we 
introduce a mass countcrtcrm Sm to cancel this divergence. Numerically it is canceled non- 
pcrturbativcly in the ratio between the three-point function and the two-point functions 
to obtain a matrix element, or in the difference between binding energies of heavy-light 
mesons. 

The bare heavy propagator on the lattice is 



a ^i-^ f —a[Sm + T^ip)] 



n 

Z2hS''{p). (4) 



By choosing the renormalization conditions 

{s^)~^{p)\iP4-^o = m, (^"i = -s(p4 = 0), 

the constant Z2h is then: 



ip4— >0 



d{ip4) 

The bare vertex function V^^ip) is defined as: 

V^ip) = (5^)-^(p)J]e^^-(^-^)(/.^(a;)Og(0)F(l/))(5^)-i(p) 

= J]e^^'-(-^)(/,«(a;)Oj(0)/i«(y))(5«)-i(p), (5) 

where 

Og(0) = Z^O5(0). 
We will see below that Vjf (p) can be written as 

v^ip) = {i + sv)uipHr'pMp) 

^ {l + SV)V,fip). (6) 

SV is given by all the IPI one-loop diagrams containing the vertex. 

We obtain {H**\0^\H) = Z-\H**\Ofj\H) where Zv = Z2h{l + SV) and {H**\Ofj\H) 
was computed on the lattice. 

This paper is organized as follows: in Sec. 2 we recall the action we use for the heavy 
quark, we give the corresponding Feynman rules and we clarify our notations; in Sec. 3 
we give the analytical expression for the heavy quark self-energy, in Sec. 4 we give the 
analytical expression for radiative corrections to the derivative operator in lattice HQET. 
We briefly conclude in Sec. 5. 
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2 Heavy quark action and Feynman rules 

The lattice HQET action for the static heavy quark is 



^HQET = ^3 ^ |/^t(^) ^h{n) - [/i'^^P(n - 4)h{n - 4)] + a6mh^ {n)h{n)^ , 



where U^^^{n) is a hnk built from an hypercubic blocking. 

We will use in the rest of the paper the following notations taken from [6]- [8]: 

Irl:^^ 



h{n) = / e*f"/i(p), 
Jp 



«'50 .a, 



2! 



J p J p 



= cos 



Fx — sin akx, 



= cos 



k, 



S/j, = sm 



iV^ = sin 



A 

In the Fourier space, the action is given at the order of C((7q) by: 

^HQET ^ f a-i/it(p)(i _ e-'P^'')h{p) + 5mh''ip)h{p) 
Jp 



+ 



2! 



5{q +p' - p){h){p)Bl{q)T''h{p')e-'^'+P'^'^-^ 

S{q + r + p'- p)h^{p)Bl{q)B\{r)T''T^h{p')e-'^P^+P'^^ 



I p J p' J q J r 

The block gauge fields -B^ can be expressed in terms of the usual gauge fields: 



n=l 

where i?^"^-* contains n factors of A. At NLO, it was shown that we only need B^^^^ [9]: 

V 

(ii = (2/3)«i(l + 02(1 + tts)), (^2 = (4/3)aiQ;2(l + 2q;3), 0^3 = 8Q;iQ;2tt3- 

Two sets of ttj's have been chosen: 1) ai = 0.75, ^2 = 0.6, a-^ = 0.3 (which has been chosen 
in our simulation and has been motivated in [2]) and 2) ai — 1.0, 0:2 = 1.0, q;3 = 0.5, 
motivated in [10]. We will label these two sets respectively by HYPl and HYP2. 

The Feynman rules can be easily deduced (they must be completed by the application 
of /i^^): 

heavy quark propagator a(l — e^*^*" + e)""*" 

vertex i;%,,(p,p') -igoT^5,,e-'^^+P'^n 
vertex V^^^f,^ggip,p') -lagi6,,{T\ T''}e-(^'4+f4)f 

gluon propagator in the Feynman gauge a^(5^,^(5"''(2iy + a^A^)~^ 

Note that p' and p are the in-going and the out-going fermion momenta, respectively. 
We also introduce an infrared regulator A for the gluon propagator. We symmetrize the 
vertex V^^,^f^gg by introducing the anti-commutator of the SU{3) generators, normahzed 
by a factor |. The gluon propagator and the vertices are defined with the A field. The 
coefficient X]i=i ^li = -^(^4) will enter as a global multiplicative factor in all the integrals 
expressed below. We have chosen the Feynman gauge: since one calculates the rcnormal- 
ization of a gauge-invariant operator, the renormalization factor Zjj is gauge invariant. 

3 Heavy quark self-energy 

From (4) we have E(p) = — (Fi -|- F2), where Fi and F2 correspond to the diagrams shown 
in Fig. 1(a) and (b): 
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4: 2 f H{Ni) e-i{k4+2ap^) 

^ ~ ~3^^° 2W + a2A2 1 - e-^C^^+ap*) + e 

1 2 / ^(^^4) e-i(fc4+2ap4) 

~3a^° A {N4 + iE){N4 -iE)l- e-^C^^+'^P^) + g 
Note that Latin indices are spatial and 



\ i i<j 

E' = 2argsh(E). 

In (9) we have ehminated properly the non covariant pole = —ap^ + ie by closing 
the integration contour in the complex plane '^{k^) < where there is the single pole 
N4 = —iE. Furthermore the integrals along the lines k^ = ivr + ik'^ are equal, because the 
integrand is 2'7r-periodic. 

Finally we have in the limit ap^ 0: 




3a 
4 



H{-iE) 



AEVTTE^ 1 - e^' 



2. f H{-iE) 



1 

+ - 



(10) 



P p+k P 
(a): Sunset diagram 




p p 
(b): Tadpole diagram 



Figure 1: : Self-energy corrections 
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We find: 



127r2 

The tadpole diagram F2 is 



Fi ^ {/i(a,)/a + zp4[2 Inla^A^) + h{a,)]] . (11) 



=ap4->0 -- — (l/a-Zp4) 



2 3a Jk 1W 

^ -Yf|i(V«-^P4)/3K). (12) 

The factor 1/2 is introduced to compensate the over-counting of the factor 2 in the Feynman 
rule of the 2-gluon vertex when a closed gluonic loop is computed. We can point out that 
the divergent part 

'^oioii) = T^fV <^o(Q;j), (7o = /i + /s (13) 

of the self-energy is smaller with the sets HYPl and HYP2 of the aj's than with the cor- 
responding contribution without "hypercubic" links [5], as shown in the Table 1: ctq (a^ = 
0) = 19.95, ao (HYPl) = 5.76 and ctq (HYP2) = 4.20, in good agreement with computa- 
tions made by the ALPHA collaboration [10], which compares the pseudoscalar heavy-light 
meson effective energy with different static heavy quark actions, and by Hasenfratz et al 
[11]. Qualitatively, one expects that the hypercubic blocking reduces the radiative correc- 
tions since it amounts roughly to introduce an additional cut-off in the integrals. 
The wave function renormalization is 

Z2H{a,) = l + ^[-2ln{a'X') + Z2{a,)], ^2 = /s - /a- (14) 

\z2 \ is also reduced by the hypercubic blocking, as indicated on Table 1: Z2 (ctj = 0) = 24.48 
[5], Z2 (HYPl) = 2.52 and Z2 (HYP2) = -3.62. 

4 Derivative operator in lattice HQET 

We have to renormalize the operator Ofj = h^jij^Djh^ . Following [6], 



«'E^5W = a'^E^''N^^^'^^N^''(^ + ^V/^''W7.7't/i(n-J>^(n-j) 

n n 

= / f a-'5{p-p')h''{p){i%lh,)h''{p') 

+ igoj I 1 5{q+p' -p)h''{p)^a''c^A%q)n^{p') 

J p J n' J n 



,2 



jJ J/(<i + ^+P'-Pyh''(P)T''T'^'^r'SjA]{q)A%^^^^^ (15) 



Note that we have chosen to not submit the covariant derivative to the hypercubic 
blocking. 

The vertex function is obtained by writing = V^j + V^j + V^j, corresponding to 
the diagrams (a), (b) and (c) in Fig. 2; V^j{ai) = u{p)'ji'j^u(j))Vj'{ai), k = 0, 1, 2. The 
contribution V^, is then given by computing 



2W + a2A2 ^'"^ ^ + ""^^^ (1 _ e-^(A^4+aP4) + e)2 



3a 



4^ 2 [ JHN^ 
3a^° A 2Vr + a2A2 

3a^° A 21^ + a2A2 



I^Fj + apj cos fcj) e~ 



e * 2 



I^Fj + aj9j cos fcj) (1 — 2ap4 



X — g-*(fc4+ap4) _|_ g 

1 



2^sin(teE±) +e 



. fc4 + ap4 



1 2 • 



(16) 



We can get rid of the integrand proportional to F^, because it is an odd term. It remains 

_ _4. 2 f HjNi) cos kj 

The integrand has poles at A^4 = ±zi?, k^ = 2zargth (|). Once again we close the integra- 
tion contour around the single pole A^4 = —iE. Since j is spatial, the "sail" diagram drawn 
on Fig. 2(b) does not give any contribution, thus V^] = 0. There is finally the tadpole 
contribution V^j which is given by computing 



K'(«,; 



11- 2 



T7 1 2 23 
12^2^^ 12.23. 



(17) 



We have finally {H**\0^\H) = Z:^\ai){H**\0^^\H){ai) where 



Z2h{a,)[l + 5V{ai)l 

ff g(iV4)cOsfc, 

V(2VF + a2A2)(2iAr4 + eM4)2 ' 4PF 



3 

4^1^ 1 [ H{-iE) cos kj 



+ 



1 



9l 



3 2 Jk {2E)^VTTE^ 127r2 



12.23 



9o 



[21n(a2A2) + /. 



:i8) 





Figure 2: : Operator corrections 
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«j = 


HYPl 


HYP2 


/l 


7.72 


1.64 


-1.76 


f2 


-12,25 


1.60 


9.58 


h 


12.23 


4.12 


5.96 


h 


-12.68 


-10.32 


-8.18 




19.95 


5.76 


4.20 




24.48 


2.52 


-3.62 


Zd 


11.80 


-7.80 


-11.80 



Table 1: Numerical values of the parameters /i, /2, /s, /4, ctq, Zd defined in equations 
(11), (12), (18), (13), (14) and (19), respectively. 

Zv{oii) = 1 + -^za{ai), za = Z2 + U (19) 

The numerical values of Zd are indicated in Table 1. 

Remark that infrared divergences appearing in and 1 -|- 6V cancel and there is 
no dependence on a, a further consequence of the /i independence of the matrix element 
{HQ{v')\h{v')^i'y5Djh{v)\H{v)) at zero recoil. Note also that as already mentioned the 
tadpole diagram of the operator vertex is not smoothed by the hypcrcubic blocking; there- 
fore it is quite large. On the other hand the tadpole contribution to the self-energy is 
smoothed by the blocking. The final result is a large positive correction to the renor- 
mahzed matrix element. By fixing qq = 1 {/3 = 6.0), this gives Z^^(HYPl) = 1.07 and 
Z-^(HYP2) = 1.10, thus the discrepancy between Z-^(HYPl) and Z-^(HYP2) is small 
(3%); without hypercubic blocking one would obtain Z]j\a, = 0) = 0.90. We can think of 
applying a boosting procedure; in the pure HYP case the boosting plaquctte correction is 
very small [8] (equation 19 and below). On the other hand we have to take into account 
that the covariant derivative operator involves links without hypercubic blocking, therefore 
one should employ a different prescription for the diagram drawn on Fig. 2(c), possibly 
leading to a larger tadpole contribution from the operator to Z^^{ai ^ 0), and therefore 
a larger positive correction. Anyhow this kind of recipe would not lead to Z^^{ai ^ 0) 
significantly larger than 1.1. 

With our exploratory lattice study and taking account Zp-'^(HYPl), we find ri(l) = 

0.41(5)(?), r3(l) = 0.57(10)(?) and tI{\) - t?(1) = 0.15(10), where systematics are un- 

2 2 2 

known; one is then not too far (within Icr) from saturating by ground states the Uraltsev 
sum rule [12] V |t^"^(1)P - |t["^(1)P = \. However the relation - nl > 9 A2t?(1) 

2 2 2 

[13], with A = Mh* - Mh ^ 0.4 GeV [1], fi^ = 0.35 GeV^ leads to ^1 larger than 0.6 
GeV^, which is significantly above experimental determination by moments. 

5 Conclusion 

In this paper we have calculated the radiative corrections to the covariant derivative op- 
erator h'^ff^Djh in lattice HQET with an hypercubic blocking of the Wilson line defining 
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the heavy quark propagator. This determines the renormahzation of the operator which is 
used to estimate the Isgur-Wise functions between the ground state and the L = 1 exci- 
tations at zero recoiL We find that there is a global, but moderate, enhancement of ti (1) 
and T3(l) with respect to the bare quantities computed on the lattice in the case where 

2 

one introduces fat timelike links, while there is a reduction with a simple Wilson line. 
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